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Abstract 

The pion charge and scalar form factors, Fi{Q^) and Fq{Q'^), are first calcu- 
lated in different forms of relativistic quantum mechanics. This is done using the 
solution of a mass operator that contains both confinement and one-gluon-exchange 
interactions. Results of calculations, based on a one-body current, are compared 
to experiment for the first one. As it could be expected, those point-form, and 
instant and front-form ones in a parallel momentum configuration fail to reproduce 
experiment. The other results corresponding to a perpendicular momentum con- 
figuration (instant form in the Breit frame and front form with = 0) do much 
better. The comparison of charge and scalar form factors shows that the spin- 1/2 
nature of the constituents plays an important role. Taking into account that only 
the last set of results represents a reasonable basis for improving the description 
of the charge form factor, this one is then discussed with regard to the asymptotic 
QCD-power-law behavior Q~^. The contribution of two-body currents in achieving 
the right power law is considered while the scalar form factor, Fq[Q'^), is shown to 
have the right power-law behavior in any case. The low-Q^ behavior of the charge 
form factor and the pion-decay constant are also discussed. 

PACS numbers: 12.39.Ki, 13.40.Fn, 14.40.Aq 
Keywords: Relativistic quark model; Form factors; Pion 



1 Introduction 

The pion has a double interest. As a physical system, it allows one to learn about its 
hadronic structure and how QCD is realized. As a test case, it can be used to check 
methods to describe properties of strongly interacting systems. Its charge form factor. 
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which has been measured up to ^ 10 (GeV/c)^ [H El El IH El E] , has thus been the object 
of a lot of attention. Noticing that the distinction is not always relevant as soon as some 
approximations are made, it has been considered in both field-theory [7| [8l[9| [T0 | [TT | [T2 | [T3] 
and relativistic-quantum-mechanics (RQM) approaches [Tl[T5l[T6t[T7t[T8l[T9l[20|[2T|[22l 
[221 [211 [25] • Most often, results have been obtained and analyzed within the standard 
light-front formalism (g+ = 0). 

Recently, more extensive studies looking at the role of various momentum configurations 
within field-theory approaches, with the transfer momentum parallel or perpendicular 
either to the average momentum of the system or to the front-orientation, have been 
made [HI [121 [IS]- A large sensitivity to various one and two-body contributions was 
observed though the total result is well determined. Not much attention was however 
given to these results, perhaps because the momentum configuration g+ = generally 
imposes itself as the most appropriate. 

The situation in relativistic quantum mechanics is somewhat different. On the one 
hand, in absence of a close relationship to field theory, there are currently discussions as to 
which wave function or/and which current will be the best to reproduce experimental data 
[151 [201 [231 [23 • The only study that tried to incorporate a more founded wave function 
[T6l [T7] failed to account for the measurements of the pion charge form factor and, in this 
order, was forced to invoke quark form factors. These ones could account for some missing 
physics which is expected to play some role like that one underlying the vector-meson- 
dominance phenomenology. On the other hand, it appeared that the estimate of this 
observable [22l [21] could be very sensitive to the form employed in implementing relativity 
|26j . quite similarly to what was observed for spinless constituents ^27]. Moreover, an 
extension of this last work [28] showed a strong sensitivity of instant and front-form results 
to the momentum configuration, like in field-theory based approaches [HI [121 [13] . It would 
not be surprising whether a similar effect shows up in the pion case too. Altogether, there 
are therefore many reasons for motivating further studies of the pion charge form factor 
within relativistic quantum mechanics. One could add that specific predictions from QCD 
relative to the asymptotic pion charge form factor [211 [301 [31] or to the Goldstone nature 
of the pion have been hardly discussed within the RQM framework. 

When looking at the pion form factor within relativistic quantum mechanics approaches, 
further work could be motivated by the numerous differences that are expected with the 
scalar-particle case considered in previous works [321 [23 [2E] • Apart from the fact that 
the relevance for a physical system of conclusions achieved for an academic one has been 
questioned, the spin-1/2 nature of the constituents represents an important if not a major 
difference as will be seen. This can affect the solution of a mass operator, the current 
(Lorentz vector or Lorentz scalar) and, ultimately, the charge and Lorentz-scalar form 
factors as well as their asymptotic behavior. Confinement represents another difference 
characterizing hadronic physics. How these differences show up when comparing form 
factors calculated in different forms (and different momentum configurations) together 
with a single-particle current is to be determined. This is our first goal. 

Anticipating that a reasonable starting point for further improvement is only given by 
the standard instant- and front-form approaches (Breit-frame and g"*" = respectively), a 
second goal is to determine the role of two-body currents in getting the right QCD 
asymptotic behavior for the pion charge form factor within these approaches [30l [31] . 
Addressing this question in a RQM framework may look premature. Apart from the fact 
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it should be ultimately considered, we notice that previous work with scalar particles 
succeeded relatively easily to reproduce the expected power-law behavior of form factors 
in the standard instant- and front-form approaches [271 [28]. This gives us some confidence 
in considering the problem here. It involves the interaction at short distances, due to one- 
gluon exchange and, of course, the spin of the constituents. With this respect and though 
there is no known scalar probe of practical relevance, the consideration of the scalar form 
factor, beside the charge one, is especially useful in revealing features specific to their 
asymptotic behavior. The contribution of two-body currents was considered in the past 
with the aim to restore the appropriate asymptotic behavior but this was for a calculation 
involving the "point form" approach and scalar constituents [33]. Moreover, the origin of 
a wrong asymptotic behavior in this last approach differs from the one considered here 
for the pion in the standard instant- and front-form approaches. 

While considering the two goals mentioned above, we found that the main trends evi- 
denced by the results were rather insensitive to various ingredients entering the mass ope- 
rator. Quantitatively, some sensitivity was observed however. We will therefore present 
here the main lines, leaving for future work a more complete discussion of the quantitative 
aspects. With this last respect, we will only mention the points that could be questio- 
nable. Preliminary results were presented in ref. [51] while the role of different forms was 
recently discussed by He et al. [21], using phenomenological wave functions. 

The plan of the paper is as follows. In the second section, we consider the description 
of the pion in different forms. We, in particular, discuss the mass operator which contains 
both a confining part and a Coulomb one. The expression of the charge and scalar form 
factors in different forms is given in the third section. One-body and some two-body 
currents are considered. The expression of the pion-decay constant, /^r, is also discussed 
as well as quantities that involve this constant, namely the charge form factor in the 
asymptotic domain and the charge radius (with some approximation) . The fourth section 
is devoted to the presentation of the results with the one-body current and to their 
discussion, for the low- as well as the high-Q^ range. The analysis of the asymptotic 
form factor is considered in the fifth section. The conclusion and further discussion is 
given in the sixth section. Three appendices contain details about accounting for the 
spin-1/2 nature of the constituents in the mass operator, the behavior of the solution 
at short distances for a strong Coulomb-type interaction, and about two-body currents 
contributing to the charge form factor. 

2 Pion description in different forms 

Prior to any calculation of form factors in relativistic quantum mechanics, two ingredients 
have to be specified. On the one hand, a mass operator, whose solutions could be used 
in different forms, is needed. On the other hand, the relation of the momenta of the 
constituent particles to the total momentum, which depends on the form, is required. 
They are successively discussed in this section. 
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2.1 Mass operator 



We here follow a previous work based on a quadratic mass operator together with a single- 
meson-exchange type interaction [27j. This approach, specialized to scalar particles, was 
able to provide a good account of form factors corresponding to an "exact" theoretical 
model in both the instant- and front-form approaches. In the case of an infinite-mass 
exchange interaction, the "exact" form factor could be reproduced, ensuring the correct- 
ness of a minimal number of ingredients for both the mass operator and the currents. In 
the pion case, two further ingredients have to be considered: the spin-1/2 nature of the 
constituents and confinement. A mass operator that accounts for these features is the 
following: 

Ml Mk) = ((2 ek + astr- Vof + A) 0o(^) 



(2^)' {k-k'f ^ 

where = y^m^ + k"^. The spinor part of the pion wave function factors out and has 
therefore been omitted in writing the above equation (see appendix lAl for some detail). 
The only but important effect due to the spin-1/2 nature of the constituents, in compa- 
rison with the scalar-particle case, is the replacement in the last term of a factor m? by 
Cfc Cfc/ {2 — rn?q {e\ e^,)/(2e^ ^\'))- The equation being written in momentum space, 
the distance r should be considered as an operator. In order to simplify the writing of 
some expressions later on, we also introduce a wave function defined differently, (pik) = 
-J&k 4>o{k). 

The confinement term together with the kinetic-energy one appears in a quadratic form, 
in accordance with the quadratic character of the mass operator. This model can 
roughly reproduce the Regge trajectories and the first radial excitations of mesons with 
a string tension ast = 0.2 GeV^ and Vd = 0.486 GeV. The other part of the interaction is 
inspired from the single-gluon exchange one. It incorporates the standard Coulomb part 
but it also involves non-relativistic corrections that are usually omitted (see appendix IM . 
It contains the QCD coupling, here replaced by an effective one, gl^^/A, and a factor 16/3 
representing the value taken by the color matrices for a meson. In accordance with the way 
we derive this part of the interaction, there is no specific contribution due to the standard 
spin-spin term. This one is actually cancelled in the present case by a term with the same 
spatial structure, generally omitted. The above effective one-gluon exchange interaction 
could a priori have a more complicated expression. In absence of a detailed study, we 
assumed the simplest choice compatible with reproducing the one-gluon exchange in the 
perturbative regime, the only parameter being the effective coupling, gl^f/A. Once the 
above ingredients are fixed, the pion mass can be reproduced by choosing appropriately 
the quantity, A. The equation can be solved in various ways. Working here in momentum 
space, we consider the operators agt r and ((Xst as the limit for G going to infinity of the 
quantities, G (1 — exp{—ast r/G)) and G^ (1 — exp{—{ast r/G)^)), and Fourier transform 
them. The solution is expected to have the high-momentum behavior, (j)o{k) oc up 
to log corrections. For a strong enough coupling, these corrections can lead to a change 
of the power-law behavior. One thus expect 4>o{k) oc k~^ for (8/3) gl^j/^ATi) = 0.5 (see 
appendix [B]). A critical regime could be reached when this quantity takes the value 2/7r 
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which is in the range of expectations when the current QCD couphng, gljf/ (4 vr) = a^, 
is used. The power-law behavior of 0o(^) is quite important as it determines that one 
for form factors [37], unless some specific suppression occurs in relation with a particular 
probe. 

For practical purposes, we use = 0.2 GeV, M^, = 0.14 GeV, = 1 GeV/fm= 
0.2 GeV^ and Vc = 0.5 GeV. For the strong QCD coupling, we distinguish two regimes, 
a low- and a high-energy one. In the first case, the effective coupling is taken as 
(8/3) g^jj-/ (An) = 0.5, which corresponds to the case mentioned above where the wave 
function should scale like at large k. This coupling value is somewhat a compro- 
mise between larger and smaller values expected respectively in the low- and high-energy 
regime (see also below for other effects). In the second case, anticipating on the fact 
that the coupling should go to 0, it was assumed that the solution behaves like A;~^/^ for 
k >5.6 GeV. As a result of these choices, the behavior of our solution slowly changes from 
the behavior k~^ obtained around 1 GeV to the one ascribed beyond 5.6 GeV. One then 
obtains the value of the last parameter, A = —0.487 GeV^. 

Being interested here in gross features, we only considered the most important parts of 
the interaction: the confinement and the one-gluon exchange, which are essential to re- 
produce respectively the dominant aspects of the hadron spectroscopy and the asymptotic 
behavior of form factors. We did not therefore try to improve upon the above model. On 
the one hand, the derivation of the effective interaction to be used in relativistic quan- 
tum mechanics is largely open. This has been done to some extent in the scalar-particle 
case where the effect of retardation was found to lead to an effective coupling 2-3 times 
smaller than the free one [38] • The estimate of this effect and other ones is likely to be 
more complicated for the one-gluon exchange interaction [SH], of interest here. On the 
other hand, the running character of the strong QCD coupling, a^, should be accounted 
for. Due to the relatively enhanced weight of high momenta in the solutions of eq. (Q, 
one can anticipate that values of this coupling smaller than those currently referred to 
for the low-energy domain should be used. An important point is that, in both cases, the 
most singular part of the interaction, which determines the asymptotic behavior of form 
factors [40j and is given in eq. ([T]) by the Coulomb term, be preserved (up to log terms). 
Of course, the interaction can be improved on a phenomenological basis, by requiring for 
instance a closer fit of the parameters to the pion radial excitations at 1300 MeV and 
1800 MeV or to the pion-decay constant, f-„. 

2.2 Constituent and total momenta in different forms: relation 

In order to determine the expression of form factors, beside a solution of a mass operator 
discussed above, the relation of the constituent momenta, pi, p2, to the total momentum, 
P, and the internal variable, k, is needed. Two ingredients enter their determination. 
The first one involves the relation of the sum of the constituent 3-momenta and the total 
momentum: 



where ^'^ characterizes the approach under consideration and, especially, the symmetry 
properties of the hypersurface which physics is described on. Most often, represents 
the orientation of a hyperplane and the above equation can then be obtained, up to a 




(2) 
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phase, by integrating plane waves, exp(i (pi +p2 — P) ■ x), over the hypersurface ^ ■ x = 0. 
The second ingredient involves a Lorentz-type transformation of the constituent momenta 
which generalizes that one introduced by Bakamjian and Thomas |41j : 

Pi,2 = ±k±w +wek, ei,2 = w^ek±w-k, (3) 

+ 1 

where = y/l + w"^ {w'^ = {w^y — w"^ = 1). Combining the two equations allows one 
to write: 

w'' = h = , (4) 

2ek 2ek ^(^ . p)2 + (4 e2 _ j\^2) ^2 + ^ . p 



where w'^ depends on the approach through that one of C,^. It is noticed that, repre- 
senting an orientation, its scale is irrelevant, what can be checked on eq. ([2]) or on the 
last expression for w'^. The expressions taken by and w'^ are given in the following for 
different forms. 

• Instant-form approach 

e oc with (A°, Ao) = (1, 0) (A^ = 1). (5) 



Front-form approach 



^'^occu^, with cJ) oc (1, fi), = 0), (7) 

w = - h- — - — , w^ = —^+ " 



2 Cfc 4:ek Ep - P -h 2 Ack {Ep - P - h) 

where n has a fixed orientation. 
• Dirac inspired point-form approach 



oc , with u) oc (1, m), (m^ = 0) , (9) 

w = — + - — ^ — , w =w^ + ^ • 10 

2efc 4:ekEp-P-u 2efc Aek{Ep-P-u) 

where m, contrary to the front-form case, has no fixed orientation. 

In the absence of interaction (2efc = M), the three expressions obtained in different 
forms for w and become identical to the standard kinematical boost for free particles. 
In this limit, one recovers the expressions relevant to an earlier "point-form" approach 

^^ocA^, with A^ = P^//P2. (11) 

^=^' ^" = % (12) 

This "point-form" resembles the Dirac point form in that the dynamical or kinematical 
character of the Poincare generators is the same. However, as mentioned elusively by 
Bakamjian [13] and explicitly by Sokolov [44j, this "point-form" implies describing the 
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E, , , p E, , Pf 

Figure 1: Graphical representation of a virtual photon absorption on a pion together with 
kinematical definitions. 

physics on hyperplanes perpendicular to the velocity of the system and differs from the 
Dirac one, which involves a hyperboloid surface. The kinematical character of boosts and 
rotations in the above "point-form" stems from the invariance under Lorentz transforma- 
tions of the hyperplane defined by the condition that f ■ a; is a constant, where v is the 
4- velocity of the system {y^ oc P^). To some extent, the property is trivial as the frame 
used to describe the system changes at the same time as this one is boosted. A related 
feature is the fact that it requires a constraint on the interaction much weaker than in 
the other approaches. When deriving a mass operator, any Lorentz invariant interaction 
is sufficient while, in the other cases, the consideration of the interaction at all orders is 
required [52] . A similar statement holds for other quantities such as the charge form factor 
at = or the pion-decay constant (see below). The above "point-form" approach thus 
possesses properties that make it definitively different from the other ones. In some sense, 
it contains the effect of the boost transformation common to all approaches but without 
the constraints attached to the fact that the hypersurface which the system is described 
on is uniquely defined and independent of its velocity. As is well known, this requires 
interaction effects, represented here by the dependent term at the r.h.s. of eq. (jl]). 

3 Current and form factors 

We consider in this section the pion form factors. There are two of them, corresponding 
to the coupling to a photon charge form factor) and to a possible scalar probe 

(Fo(Q^): scalar form factor). Their definition we refer to here may be found in ref. [27] . 
Considering both of them can give a better insight on the implementation of relativity, 
especially with respect to the non-zero spin of the constituents. The studies are performed 
on the basis of a single-particle current. Two-body currents will be considered but only 
for those approaches that already provide a relatively good account of experiment. We 
therefore discard large contributions from two-body currents that could be required to 
preserve properties related to the Poincare space-time translation invariance [2HI SS] • In 
short, these last currents account in a RQM framework for the equality of the momentum 
that is separately transferred to the constituents and to the whole system, which is fulfilled 
in field-theory approaches in a straightforward way. 
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3.1 One-particle current 



The form of the single-particle current to be used with the above mass operator has been 
discussed at length in ref. [27] for spinless constituents. Interestingly, it has been found, 
since then, that form factors in different forms, due to minimal consistency requirements, 
could be written in a unique way p8]. They should be completed to account for the 
spin-1/2 nature of quarks. Referring to fig. [1] for the kinematical notations, they now 
read: 



dp ^f{Pf+p) ^i-{Pi+p) ,//P/-Pn2 



(27r)3 ep(20-P/) i2^^■Pi 



X Upf + Pf {Pi + Pf 



^1/4 



0o(( 



Pi- P. 2 



S 



{Pf+P^ + 2p)-{if + ^i) 
In this equation, and 5* result from the following trace over 7 matrices 



cfo .(13) 



S 



4Tr 
2 



75 



(mg + 7 ■ Pi) ^ (mg + 7 - pj) (mg - 7 • p) 



1 



75 



2 ' 2 

Pt {Pf ■P)+P'f {Pi ■p)-p'' {Pi ■ Pf) + ml {pf +Pf+ p^) 



2 

,2 /^i^ 



4 Tr 
2 m, 



75 



(mg + 7-pj) (nig + 'y-pf) (^g - 7 ■ p) 



75 



Pf ■ p + Pi ■ p + Pi ■ Pf + m" 



(14) 



In the case of a point-form approach inspired from the Dirac one, it turns out that 
expressions for form factors can be obtained from the above ones in the front form by 
integrating them over the orientation of the front together with an appropriate weight. 
They can be obtained from a straightforward generalization of those for the scalar-particle 
case [28] . 

It can be checked that the expression of the charge form factor at = considerably 
simplifies. Making an appropriate change of variable, it reduces in all cases to a unique 
expression: 



Fi'\Q' = 0) = J 



dk 
(2^ 



1 



(27r) 



efc 



1. 



(15) 



We stress that the result is independent of both the momentum of the system, P, and the 
orientation of the hyperplane given by a property which is not trivial in most cases. 
With this respect, the presence of the factor, {pf + pj + 2p)-^, at the denominator in 
eq. (fT3|) is essential. This one accounts in a hidden way for some two-body currents [27] . 
It is also noticed that the above expression of the charge form factor at = is fully 
consistent with the expression of the normalization of the solution, (f)o{k), of the mass 
operator, eq. ([1]), which involves the same integrand. 

In the case of the scalar form factor, and for the instant form, it is generally found that 
its value at = depends on the momentum of the system. In other words, it does 
not verify a minimal Lorentz invariance property. This one can be fulfilled by introducing 
a correction factor in the expression of the scalar form factor, FqIQ"^). Denoted cfo, 
this quantity amounts to account for two-body currents. Its expression, whose form 
is suggested by the consideration of a simple interaction model, has been considered 
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independently [25] . Symmetrizing its effect between initial and final states, it is taken as: 



c/o = 1 + 2 ioiki) + g{kf)) e. 



p 



(e/ + ep){ei + Cp) - Ef Ej 
(e/ + ep)(ej + ep)(e/ + ei)' 



(16) 



where the function g{k) can be obtained from the following quadrature: 



e, {4el-M') <Pl{k) = / 

Jk 



■oo 



dk' k' Cfc' 00 (^') • 



(17) 
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In the case of a scalar-particle model together with a zero-range interaction, g{k) is equal 
to 1 and the above factor c/o then allows one to reproduce the "exact" scalar form factor. 



A few remarks are in order here. A first one concerns the Lorentz-invariance properties 
of form factors in different forms. While the point-form approaches evidence this property, 
the instant- and front-forms ones do not. These last approaches can thus be considered 
for different momentum configurations which include the standard ones (Breit frame and 
g+ = respectively) but also non-standard momentum configurations like a parallel one 
where the transfer momentum is oriented along the average momentum of the system 
((Pj — Pf) II {Pi + Pf) II n; Ei 7^ Ef). Their comparison could be instructive but this 
supposes that no other symmetry is significantly broken at the same time, which has to be 
checked. A second remark concerns the expressions of form factors in the front-form case. 
The current ones, generally given in terms of the Bjorken variable, x, and the transverse 
momentum, are recovered from eq. ( fT3ll by making a change of variable. The last remark 
concerns the calculation of form factors in the earlier "point-form". As the initial and 
final states have generally different momenta, two different 4- vectors, = P-^/M and 
= Pf/M, and therefore two different hyperp lanes, are then involved [17]. 

When comparing the various approaches, an important feature emerges. The boost 
transformation allowing one to get the wave functions describing the initial or final states 
from the solution of a mass operator only depends on the constituent mass in the standard 
instant- and front-form approaches [471 EH] while it also involves the mass of the system 
in all other cases. In a few of them considered later on (point-form and non-standard 
approaches in a parallel momentum configuration together with \Pi + Pf\ —>■ oo), form 
factors in the single-particle current approximation are shown to depend on the momen- 
tum transfer Q through the ratio Q/M [28]. This has striking consequences in the case 
where the mass of the system is small in comparison with the sum of the constituent ones, 
like for the pion. 

3.2 The pion-decay constant, f.j^ 

At first sight, the calculation of the pion form factors can be performed independently of 
any knowledge about the pion-decay constant, fj^ (~ 93 MeV in our conventions). This 
quantity is nevertheless relevant here as it enters in the QCD prediction of its asymptotic 
behavior for the charge one. It is successively considered in the instant and front forms. 

• Expression of f^^ in the instant form 

In a first approximation, the pion-decay constant in the instant form could be obtained 



FoiQ'). 
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from considering the matrix element of the axial current between the pion and the vacuum: 

f /m„ ~ N ei + 62 



xTr 



(27r)3y V " 2ei 62 

(mg + 7 ■ pi) 7/^ 75 (mg - 7 • P2) 75 



where the expression of k is given at this point by k'^ = (pi + j92)^/4 — m"^. Apart from 
the fact that the consideration of time and spatial components give different answers, 
which is not a surprise in a non-covariant approach (this is reminded by a question mark 
in front of the r.h.s. of the above equation), it is noticed that, in the cm., the l.h.s. is 
proportional to the pion mass while the r.h.s. is not. A similar problem is encountered 
for the charge form factor, revealing striking features for a strongly bound system [27]. It 
points to a missing contribution in the time component of the current given by the trace 
at the r.h.s. of eq. ( fT8l) . This one, given by (ci + 62)/^^, should be completed by the 
quantity [E^^ — ci — e^j/mq, which can be seen to be an interaction term and provides a 
two-body current using eq. ([1]). The expression of /^r so obtained is independent of the 
component of the current and is given by: 

By making a change of variable given by eqs. ([31 [6]) , it is found to be equal to: 

which is independent of the momentum of the pion, P, and therefore Lorentz invariant. 
• Expression of /^r in the front form 

The expression of /tt in the front form in terms of the Bjorken variable and the transverse 
momentum can be found in different works. With our convention for the wave function, 
4>o{k), it reads: 



y/S /■ ,2, dx 



[^T^Y J Jxil- x) Jml + k^, 



q 



where the argument of the wave function, k, is defined by k'^ = (mg + A;^)/(4x(l — x)) — m^. 
Not surprisingly, the above result can be recovered from the instant-form expression, eq. 
([18]), in the limit P 00 and, of course, from eq. (l20l) by making the extra change of 
variable, k^ = (1 — 2x)^ {rn^ + k\)/{Ax{l — x)). One has therefore: 

h = f'/ = f^. (22) 
• Expression of for an arbitrary hyperplane 

An expression of fj^, which is generalized to an arbitrary hyperplane of orientation ^'^ and 
evidences the main ingredients, is given by: 



X Tr 



(mg + 7 ■ Pi) ^ ■ 775 (mg - 7 ■ ^2)75 



4mg 



^ ■ (Pl + P2) ^ Cfc 



m) ■ (23) 
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Pf 



Figure 2: Graphical description of a virtual-photon absorption on a pion with evidencing 
the contribution of one-gluon exchange representing the Born amplitude: kinematical 
definitions. A similar diagram with the gluon exchange in the initial state should be 
considered. 

It simplifies to read: 



The structure of this last expression is somewhat similar to that one used for the charge 
form factor, eq. f[T^ . at = 0. Like for this one, it can be checked that the result 
evidences the property to be both independent of the velocity of the system (therefore 
Lorentz invariant) and of the orientation of the hyperplane, By performing a change 
of variable, it can be cast into the form of either eq. (120|) or eq. (!2T|) . 

3.3 Two-body currents 

We here consider two-body currents that are required to recover the full Born amplitude 
depicted in fig. [21 taking into account that part of this diagram with a positive-energy 
quark is already included in the contribution of the single-particle current previously 
discussed, eq. f|T3l) . This will be done for the pion charge form factor where the problem 
of recovering the QCD prediction for the asymptotic behavior [30l |3T] is the most crucial. 
This one is given by: 



The problem is examined in both the standard instant and front forms, which are the only 
ones to provide a reasonable starting point for considering the contribution of interest in 
the present work. The expression of the two-body currents is given here in the simplest 
cases while the derivation is considered in appendix O Two-body currents could also be 
considered for the scalar form factor but, as this one turns out to have the right asymptotic 
power-law behavior, the necessity for studying them is not as strong as for the charge form 
factor. 

• Expression of F^\Q'^) in the instant form (standard Breit frame: Pf = —Pi = Q/2) 
The two-body contribution to the charge pion form factor in the instant form may be 




(24) 




(25) 
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expressed as a double integral over the momenta of that quark in the initial or final states 
which does not interact directly with the external probe (see fig. [2]). It reads: 



X 



dpi2 eii + ei2 /TUq - 

To o I i 

1 2 ("^PS ~'~ ^1 



dpf2 e/i+e/2 /mg ~ 
(2 7r)3 2e/i e/2 ^e^ 



4:71 tts I 



/i^ - iPi2-Pf2) 



pv 



(term : i ^ f) 



(26) 



where quantities that could be related to the pion-decay constant, /^r, have been purposely 
factorized. The gluon propagator is made apparent. To keep track of its effect in formal 
developments, we ascribe to the gluon a mass, /i, which is set to zero is actual calculations. 

The various quantities at the last line refer to an intermediate particle with negative 
and positive values of the quantity and to the pseudo-scalar and pseudo-vector parts 
of the Fierz-transformed gluon-exchange interaction. They are respectively represented 
by superscripts — and + and subscripts ps and pv. As for the one-body part, the matrix 
element of the current is divided by the quantity representing the sum of the kinetic 
energies of the constituents in the initial and final states, en + e/i + 6^2 + e/2- They also 
contain factors relative to the intermediate quark propagator. The summation over the 
quark spins has been performed. Their expressions read: 



1 



ps 



x2 



en + e/i + ei2 + e/2) 2 (e) + en) 

Pii-Pi2 Pj^ -Pii-Pj P%+Pi2-Pf Pn+^l iPj^+Pa+Pa) 



ml+Pfi-Pf2 



X 



{en + e/i + ei2 + e/2) 2 ej {ej + en) 
Pf-{Pfi+Pf2) {Pn+Pa) + {Pii+Pi2)-{Pfi+Pf2) Pj^ -pJ-{Pii+Pi2) (P/l+P/2 



+Pi2-(P/l+P/2) Pn -Pil-{Pfl+Pf2) Pa + ("^g +PirPi2) (P/l+P/2 



(ea + e/i + ei2 + e/2) 2e} (eii + e} + 2ei2) (e} + e,2 + e/i + e/2) 



X 



m^+P/rP/2 (m^+pi2-P/)) + (^g+P/i-P/2) 



m^+pj-pi2 J 
X {pn-Pi2 pf^ - Pii-P'} Pa + Pi2-P~f Pn + ipf^+P^ti+Pa] 



(27) 



with e 



/ 



ml+pf, pf^ 



±e 



/, Pf =P 



The above appearance of contributions referring to the negative-energy component of the 
intermediate particle in fig. [2], 1'^ and 7^^, is evident as they involve degrees of freedom 
that are not explicitly accounted for in a RQM approach. From their derivation, it can 
be checked that they vanish at zero momentum transfer. The separation into pseudo- 
scalar and pseudo-vector contributions is justified by the different role they play as for 
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the asymptotic behavior of the pion charge form factor and, to some extent, for the charge 
radius. 

The presence of contributions referring to the positive-energy component of the inter- 
mediate particle is less trivial. As the spin part of the pion wave function in the RQM 
formalism assumes a unique form, oc « 75 ^(^2)5 it cannot fully account for that part 
of the Fierz-transformed interaction involving pseudo- vector currents. The discarded term 
has an off-energy shell character however and, thus, can contribute to the form factor. 
The corresponding contribution involves both pseudo-scalar and pseudo-vector terms, J^^ 
and J^. These ones do not vanish at zero momentum transfer but their sum does as a 
result of a close relationship. For this reason, only their sum has been given in the above 
expression of two-body currents. At non-zero momentum transfers, it turns out that 
the cancellation still holds but this result is specific to the Breit frame considered here 
{Pf = —Pi = Q/2). An expression valid for an arbitrary frame is given in the appendix 

o 

We stress that, apart from neglecting higher-order terms in the coupling, a^, the above 
two-body contribution is entirely determined by recovering the Born amplitude shown in 
fig. [2] after it is added to the single-particle contribution. Getting this amplitude with 
the right strength therefore supposes that the contribution is correctly calculated as far 
as some minimal ingredients are concerned. 

From a rough examination of the above equations, one can convince oneself that the 
large limit of Fi^\Q'^) in the Breit frame will take the expected form given by eq. 

except perhaps for some factor. This one requires some care, especially with the 
treatment of the components of the constituent momenta perpendicular to the momentum 
transfer, Q. Contrary to a naive expectation, these ones can be shifted from their zero 
value by an amount approximately given by Q (/c/cfc). In order to get a correct expression 
of the asymptotic limit, it is appropriate to make a change of variable: pi, P2 —>■ k, P. 
One then obtains: 

X 



(2 7r)3 J Ck, ^Ck, 

\/3 r dkf /Trig ~, ,\ 4 



(27r)3y ek, ^ekj > 9 (i_g.iWi + g.^ 



We notice that this expression contains factors depending on the angle of Q and k, which 
apparently have no counterpart in the standard expectation, eq. ( l25i) . As this one was de- 
rived in a front-form approach, we will come back to it after considering the corresponding 
expression of Fi'\Q'^). 

• Expression of F^\Q'^) in the standard front form (q^ = 

The expression of the two-body contribution in the standard front-form approach (g+ = 0) 
is expressed in terms of the variables commonly used in this formalism, the Bjorken 
variables x and the transverse momenta, k±. It reads: 

d^ki\ dxi /rria ,\ 1 /" d'^kn dxi 

An a. 



^ {ips + + Its + It) + (te™ : ^ /) , (29) 
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where the quantities corresponding to the gluon propagator, and to the current, /^^ + 
Ipv + Ips + Ipv^ are given by: 



+ . . . = ^2 _ (p.2_p^2)^ 

u'^ + ml — i 1 + 2 (xi kf—Xf ki) ■Qj_ + XiXfQ 



-^PS ^pV ^ 5 



, Xi (xi ki±-Q± 
I + I = 

ps ' pv , / 

2[ml + [ki-Xi Q 



(30) 



Two results stem from the above two equations. The contribution of the two-body current 

to the charge form factor vanishes at = ^s expected from the fact that two-body 

currents should not contribute to this quantity in a RQM approach. The expression of 

the asymptotic form factor can be obtained by assuming that the wave function is peaked 

at low values of the k variable, allowing one to discard k terms in the matrix element 

of the current and the gluon propagator. One thus successively gets -|- ^ 1/2, 
/i^ -)-... ~ Xj x/ ancl; 

X 



^ / J^h^ fn^i ^(fc,)) 1 (31) 

{271^ J 2xf (l-Xf) ^ekf^ ^'^ 9 XiXf ^ ' 
• Comments 

The asymptotic expressions in the standard instant and front forms, eqs. (!28|) and (I3T1) . 
differ from each other. By making an appropriate change of variable however, it can be 
checked that they coincide. On the other hand, the presence of the factors 9/4(1 — Q • 

-^)~^ (1 + Q ■ —)~^ in eq. fl28l) and {9xi Xf)~^ in eq. fl3T]) makes them apparently 

different from the usual asymptotic expression, eq. (1251) . Looking for an explanation 
of this possible discrepancy, it was found that the above expectation assumes that the 
integrand entering the expression of the pion-decay constant varies like x{l — x). In 
such a case, which supposes that our wave function (/>o(/c) exactly scales like e^^''^, a full 
agreement is recovered. Actually, our expressions are more complete with some respects. 
For a part, the wave function contains some log corrections due to higher-order effects in 
the interaction. As for the extra factors mentioned above, they were obtained in the past 
on a different basis (see for instance ref. 



The two-body contribution to the pion charge form factor assumes a quite different 
expression in the standard instant- and front-form approaches. Similarly to the one-body 
contribution, where we could find a common expression depending on the orientation of 
the hyperplane which the physics is formulated on, a common expression can be obtained 
for the two-body part. This one, which can be used for further investigations, is given in 
the appendix O With this respect, it is noticed that the asymptotic behavior is always 
produced by the pseudo-vector pseudo-vector term of the Fierz-transformed one-gluon 
exchange interaction, somewhat similarly to what has been shown for a long time in 
field-theory based approaches [8], but, depending on the RQM approach, it originates 
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from that part involving the negative-energy component of the intermediate particle with 
momentum p*'^ (instant form in the Breit frame) or from the positive-energy one (front 
form with g"*" = 0). Interestingly, in an instant-form calculation away from the Breit frame 
but preserving the condition Ei = Ef, or equivalently Q -L {Pi + Pf), the part involving 
the negative-energy component tends to vanish with an increasing average momentum 
{\Pi + Pf\ —>■ oo). The asymptotic contribution is then obtained from the other part 
with a positive-energy component. This result confirms the observation that a front-form 
calculation is closely related to an instant-form one in the above large momentum limit. 

• Expression of the charge radius 

An expression of the squared charge radius is currently referred to in the literature: 

= = 0-342 fm^ (32) 

It has been derived in different ways [i9l [50l [51] . In the original work by Tarrach for 
instance, it was obtained from the pion-quark-antiquark amplitude determined by the 
coupling constant, (yf^gg. The relation to the above result assumes the equality l/fn = 
g-irgq/^q- Being sometlmes presented as a consequence of the Goldstone-boson nature of 
the pion, the question arises of whether this result could be recovered in a RQM approach. 

Examining the derivation of the above approach in an instant-form approach (Breit 
frame), we found that a half is contributed by what corresponds to the Darwin- Foldy 
contribution to the form factor in the non-relativistic limit (— Q^/Sm^), which involves 
positive-energy spinors, and the other half by a similar term involving negative-energy 
spinors accounted for by two-body currents. There are other contributions which respec- 
tively increase and decrease the above ones. They however cancel each other and have 
therefore no effect on the total result. These various results could thus be usefully com- 
pared to those obtained in RQM approaches by isolating the appropriate contributions. A 
possible problem may concern the overall factor in eq. ( 132|) . which does not appear 
explicitly in these approaches. It has to be hoped that the description we are using for 
the pion is numerically close to that one determined by the coupling g^^qg. 



4 Results in the one-particle current approximation 

We present in this section results for both the charge and scalar pion form factors cal- 
culated in the single-particle current approximation and for different RQM approaches. 
These ones include the instant and front forms with a standard momentum configuration, 
Breit frame in the former case (denoted I.F. (Breit frame)) and g"*" = in the latter one 
(denoted F.F. (perp.)). Other frames are considered with a parallel momentum configu- 
ration (denoted I.F. -|- F.F. (parallel)). Results also include an approach inspired from 
Dirac's point form (denoted D.P.F.) as well as an earlier "point form" that has been em- 
ployed in previous works (denoted "P.F."). Contrary to the other results, these ones are 
frame independent. As mentioned in sect. 3, expressions of form factors can be obtained 
by generalyzing those given elsewhere for the scalar-constituent case [28| H6]. 

Numerical results are shown altogether in figs. [3] and H] for the charge and scalar form 
factors respectively. Two ranges of momentum transfers are considered in each case, low 
g2 on the l.h.s. (0 - 0.2 (GeV/c)^) and high on the r.h.s. (0 - 10 (GeV/cf). The 
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Figure 3: Pion charge form factor, at small and high (left and right panels 

respectively). In the latter case, the form factor is multiplied by and represented on a 
logarithmic scale to emphasize the asymptotic behavior. 



first one is expected to evidence a sensitivity to the charge (or scalar) radius while the 
second one a priori emphasizes the asymptotic behavior. As this one should behave like 
Q~^, the corresponding form factors are multiplied by so that the product be close to 
a constant in the limit of large momentum transfers (up to log terms). Measurements for 
the charge form factor are also shown in fig. [31 

The comparison of different theoretical approaches evidences features quite similar to 
the scalar-constituent case [28]. The Breit-frame instant-form results slightly overshoot 
the front-form ones with g"*" = 0. The other instant- and front-form results with a parallel 
momentum configuration, which coincide with each other, show a fast drop off at small 
Q^, a property that is also found for the point-form results. As explained elsewhere [22], 
this drop off points to a squared charge radius that is determined by the inverse of the 
squared pion mass and is larger than experiment by an order of magnitude. For a part, 
this sensitivity to the pion mass explains the large discrepancies between results shown 
at the r.h.s. of figs. [3]and|ll which roughly scale like in the scalar-constituent case. 

The comparison with experiment, which is only possible for the charge form factor, 
shows that the standard instant- and front-form approaches provide reasonable results as 
a starting point while all the other ones are off, sometimes by orders of magnitude. This 
is again somewhat similar to the scalar-constituent case where the role of the experiment 
is played by an exact calculation. As has been suggested [2B] and shown later on [36|, the 
largest discrepancies point to a violation of Poincare space-time translation invariance, 
which requires specific two-body currents to be restored. This does not necessary imply 
that the other results in the standard instant- and front-form approaches are fully under 
control. 

Looking carefully at the corresponding product of the charge form factor with Q^, 
shown at the r.h.s. of fig. [3] , it is found that the appearance of a plateau, which could 
suggest that the expected asymptotic behavior is reached, is misleading. The comparison 
with the scalar form factor, shown at the r.h.s. of fig. HI is here very instructive. In 
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Figure 4: Same as for fig. [3] but for the scalar form factor Fo[Q'^). 



this case, the same product increases, in agreement with possible log(Q^) deviations and, 
moreover, is significantly larger. Globally, it thus appears that the above charge form 
factors at high Q"^ are suppressed with respect to the scalar ones, which can be better 
seen by examining results for much higher momentum transfers (see next section). The 
origin of the suppression, roughly given by a factor Q^, resides for a part in the coupling 
of quarks to the external probe, as can be checked from considering matrix elements 
of the current between positive-energy spinors [22]. A similar suppression, but within 
a truncated field-theory light-front calculation, was mentioned in ref. [52]. This result 
represents an important difference with the scalar-particle case. 

At this point, a few further remarks can be made. The spin-1/2 nature of the con- 
stituents also shows up in results at small Q^. Its effect in the better cases (standard 
instant- and front-form approaches) is found to be very similar to the one produced by 
the Melosh transformation in other works [15], with an increased contribution to the 
charge radius (Darwin- Foldy term). Results evidence some sensitivity to the dynamics 
like the description of the confinement. This has not been considered in detail here but 
we notice that reproducing the pion-decay constant allows one to reduce the uncertainty. 



5 Results with two-body currents and asymptotic 
behavior 

We consider in this section the pion charge form factor with including the contribution 
of two-body currents. This is done for both the instant-form and front-form approaches, 
respectively in the standard Breit-frame and so-called g+ = configurations. As already 
mentioned, the other approaches suppose dominant contributions from specific two-body 
currents restoring the Poincare space-time translation invariance. These ones could be 
determined along lines developed in ref. [16] but it is not clear how they will affect the 
contribution of two-body currents of interest here (possible double counting). 

As two-body currents considered in this section are mainly motivated by the asymptotic 
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1=0.2 GeV, M=0.14 GeV m=0-2 GeV, M=0.14 GeV 
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Figure 5: Form factors -Fi(Q^) and Fo(Q^) in the ultra asymptotic domain: contribution 
from the single-particle current (s.p.a.) and, for Fi{Q'^), contribution incorporating two- 
body currents (with MEC). Calculations are performed in the standard instant and front- 
form approaches. 



behavior of the charge form factor, we first present results relevant to this feature. These 
ones, multiplied by Q^, are given in fig. [5] (l.h.s.) for a range of extending from to 
200 (GeV/c)^. For the sake of the discussion, we also present in this figure (r.h.s.) results 
pertinent to the scalar form factor. The two-body currents also contribute at low and 
intermediate values of Q^. The corresponding results for the charge form factor alone are 
shown in fig. El The range under consideration is the same as in fig. [Sj but a linear scale 
instead of a logarithmic one is adopted for the r.h.s. part. 

We first notice that, contrarily to what the consideration of the r.h.s. of fig. [3] would 
suggest, the asymptotic behavior of the charge form factor is not reached yet. Examination 
of fig. [5] (l.h.s.) clearly shows that the product of with the charge form factor calculated 
from the single-particle current, after evidencing a maximum around 10 (GeV/c)^, begins 
to slowly decrease. The difference with the scalar form factor shown in the r.h.s. of 
the same figure, which has the right asymptotic behavior, is striking. Accounting for 
the contribution of two-body currents to the charge form factor allows one to obtain its 
asymptotic behavior which is expected to be given numerically by Fi{Q'^)q2^^ = 0.17 
{GeV/cY/Q^ for our pion description. As it can be seen on fig. [5] around 200 (GeV/c)^, 
this value is not far to be reached in the instant-form case but still rather far away in 
the front-form one. The different behavior can be explained as follows. The instant- 
form two-body current implies a contribution originating from the pseudo-scalar part of 
the interaction [Ipg), which cancels the single-particle one, leaving as a neat result the 
contribution due to the pseudo-axial vector part (/p^), while this destructive contribution 
is absent in the front-form case. A few further comments may provide better insight 
on the above results. First, in spite of different formalisms, present results qualitatively 
agree with those of earlier works [H [5^ . Some of the quantitative differences, especially 
the onset of the asymptotic behavior at larger in this work, can be related for a part 
to our pion description that involves sizeable high-momentum components. Second, the 
difference between the asymptotic value of the form factor given above in the text and the 
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Figure 6: Effect of two-body currents on 
intermediate momentum transfers: results 
front-form approaclies. 
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the cliarge form factor, Fi(Q^), at low and 
are presented for the standard instant and 



current one given by eq. fl25l) . 0.11 (GeV/c)^/(5^, can be ascribed to the non-perturbative 
effects mentioned in the comments about this last result (end of sect. 3.3). 

Considering now the results for the charge form factor at very small transfers (l.h.s. of 
fig- E]), it is observed that the contribution of the two-body currents tends to remove the 
difference between the instant- and front-form form factors. Such a result is not totally 
unexpected as, for a part, the two-body currents tend to restore Lorentz invariance in 
first place and perhaps some Poincare space-time translation invariance in second place. 
It actually holds for the contribution of the pseudo-scalar part of the interaction alone, 

in eq. (I28l) . what it should as the part related to the pseudo-axial vector part, 
is irrelevant for the argument. It remains that some contribution could be due to the 
confinement interaction but its non-perturbative character does not allow one to determine 
how much it affects /~ in the instant form. 

To some extent, the above Lorentz- invariance argument applies to the form factor in 
the intermediate range, shown in the r.h.s. of fig. El but the effect is larger than needed. 
We observe that the contribution of two-body currents tends to make the instant-form 
form factor closer to experiment while it makes the front-form one further away. To better 
understand the role of two-body currents depending on the form, it is useful to consider 
separately the two contributions in the instant form, 1'^ and eq. ( l28i) . The first one, 
which has no well determined counterpart in the front-form case, is negative over the full 
range of considered in this work and decreases faster than Q~^. The second one, which 
is positive and provides the asymptotic behavior, is close to the contribution Ipg + Ipy 
calculated in the front-form approach. 

Having considered two-body currents, we can now discuss the pion charge radius to 
which they can contribute. This is done in relation with the expression = 3/(4 7r^ /^), 
which gives the value = 0.342 fm^ for = 93 MeV. A value that could be more relevant 
for the discussion is = 0.263 fm^, which corresponds to fj^ = 106 MeV obtained with 
our pion description. In the instant form, it is expected that part of the contribution 
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should have a counterpart in the piece of the two-body currents associated to the pseudo- 
scalar component, We first notice that the squared matter radius has a relatively 
small role, of the order of 0.08 fm^ in the present calculation to be compared to the 
measured value of the squared charge radius, 0.43 fm^. We also notice that the part 
of two-body currents due to the pseudo-axial vector component of the interaction, also 
small (— 0.02fm^), is irrelevant for a comparison with the above prediction which assumes 
a pseudo-scalar component and a point-like pion. As already mentioned, this approximate 
prediction should be compared to that part in our calculations which corresponds to the 
Darwin-Foldy term. Isolating this term in the instant-form results, we find 0.15 fm^ for 
the one-body current and 0.08 fm^ for the two-body one. The sum, 0.23 fm^, compares 
to the expectation 0.263 fm^. In the front form, where we only have a contribution from 
the one-body current, a similar procedure gives 0.29 fm^, which also compares to the 
expectation 0.263 fm^. Thus, there is a reasonable agreement of the present results for 
the charge radius and the approximate expectation advocated in many works. Another 
relation, suggested by the analysis of various contributions to the above expression of r^, 
is worthwhile to be mentioned. The usual Darwin-Foldy correction to the form factor in 
the instant form takes the form — (5^/(8 m^) (non-relativistic limit). When the similar 
term associated here to the two-body current is considered, one gets instead — (5^/(4 m^). 
This last term can be seen as the first term in the expansion of the vector-meson 
dominance model of form factors, —Q"^ /my, where my = Am"^ 

When considering calculations of the two-body current contribution, a large sensitivity 
to various ingredients was found. Most of it points to the role of higher-order effects in 
the interaction which, for consistency, were neglected. These effects nevertheless show 
up in looking at some of the results presented here. As an example, it is likely that the 
charge and the scalar form factors in the asymptotic regime should be comparable, up to 
a factor 1-2. Examination of results presented in fig. E] rather suggest that they differ 
by a factor 10-20. This can be largely explained by non-perturbative corrections to the 
wave function entering the calculation. Thus, these effects lead to log{Q'^/m'^) terms that 
enhance the scalar form factor in the asymptotic domain while they are ignored in the 
present derivation of two-body currents based on a perturbative approach. An other hint 
about these non-perturbative effects is provided by the comparison of the contribution 
of two-body currents in different forms. While the contributions of the single-particle 
current in the standard instant- and front-form approaches are relatively close to each 
other, those for the two-body currents significantly differ. Most of the discrepancy can be 
traced back to the term whose role exceeds the restoration of Lorentz invariance at 
soon as the momentum transfer increases. 



6 Discussion and conclusion 

We first considered in this work the single-particle contribution to the charge and pion 
form factors in different forms of RQM approaches. These ones include the standard in- 
stant and front forms (Breit frame and the momentum configuration g+ = respectively), 
the same forms with a "parallel" momentum configuration and the point form. The in- 
teraction entering the mass operator is chosen as the sum of a confining interaction and a 
gluon-exchange one, which is a priori essential for getting the right asymptotic behavior 
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of form factors. For the case of the charge form factor, where measurements are available, 
the standard instant- and front-form approaches provide results relatively close to them. 
All the other results fall far apart, as a result of their dependence on the momentum 
transfer through the ratio v = (5/(2 M^-). Thus the pattern evidenced by these results is 
very similar to that one obtained with a model of scalar particles |2H] • As has been shown, 
the last results evidence a strong violation of Poincare space-time translation invariance 
while this property is approximately fulfilled by the standard instant- and front-form ones. 
This discards the underlying approaches as efficient ones to describe the pion properties. 

While there is some similarity of the above results with earlier ones for a theoretical 
model with scalar constituents, there are nevertheless differences. The most striking one 
concerns the ratio of the charge and scalar form factors. Whatever the approach, the first 
one is suppressed compared to the other one at high Q^, preventing one from reproducing 
its expected asymptotic behavior, Q~^. This result, which seems to characterize relativis- 
tic quantum mechanics approaches with a single-particle current, implies the spin-1/2 
structure of the quarks. This conclusion is comforted by the consideration of the scalar 
form factor which, up to a finite coefficient, has the right power-law behavior. 

We then looked at the two-body currents whose contributions are necessarily required 
to reproduce the asymptotic behavior of the charge pion form factor expected from QCD. 
This has been done in both the standard instant and front forms, which are the only 
approaches to provide a reasonable starting point for further improvement. The two- 
body currents have been derived on the basis that adding their contribution to that one 
generated by the wave function should allow one to recover the full Born amplitude. While 
the description of hadron properties at low- and high-momentum transfers are often consi- 
dered as disconnected, we showed that the above two-body currents were able to reproduce 
the expected expression of the pion charge form factor in the asymptotic regime. This 
result involves the pseudo-vector pseudo-vector part of the Fierz-transformed one-gluon 
exchange interaction, in agreement with what is expected in field-theory based approaches 
[8]. Quantitatively, the relevance of this result appears at quite high momentum transfers, 
probably beyond the range where the pion charge form factor has been measured. At 
very low-momentum transfers, the contribution of these two-body currents depends on 
the form. While it increases the front-form form factor, as at high Q^, it decreases the 
instant-form one. The total effect compares to the difference of the instant- and front- 
form results in the single-particle current approximation, which it tends to cancel. Such 
a result is expected for a part as far as the two-body currents we are considering should 
contribute to restore Lorentz invariance. Between the very low- and very high-Q^ domain, 
it is difficult to attribute a specific role to the above two-body currents, but it does not 
seem to affect much the comparison with experiment. 

We have not made any attempt to provide a better description of the measured pion 
charge form factor, being mainly concerned by the determination of an approach where the 
bulk contribution would be given by a one-body current so that the derivation of two-body 
currents of interest here be simplified as much as possible. Following refs. [121 |T7], the 
introduction of some quark form factor could reduce the discrepancy with experiment. 
One can nevertheless guess that a smaller quark mass would contribute to get both a 
better value of /tt and a better charge radius, hence a better form factor at very small 
Q^. A pion description with a smaller amount of large momentum components could 
also have a similar effect. As, apart from the quark mass, the parameters entering the 
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mass operator we used are essentially fixed, the only issue would be to modify the form 
of this operator. One can expect for instance that retardation effects, ignored here, could 
contribute to reduce the strength of the one-gluon exchange at short distances [39l [38] 
and, consequently, the interaction at high momenta. Thus, while relativistic quantum 
mechanics is not the most fundamental approach for describing the pion properties, getting 
these ones relatively well in this framework does not seem to be out of reach. 
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A Accounting for the spin of the constituents in de- 
riving the mass operator 

The simplest expression of the RQM pion wave function reads; 

TTl 

ip{pi,P2) oc u{pi)'y5v{p2) — Mk), (33) 

where the variable k is related to the constituent momenta by the relation 4:e1 = (pi+P2)^- 
The spinors are normalized as follows: 

2^ u[p)u[p) = — , 2^ v{p)v{p) = — ^. (34) 

spins ^"^9 spins ^"^9 

Other expressions of the pion wave function reduce to the above one, eq. (!33|) . taking into 
account that the on-mass shell character of particles in the RQM framework allows one 
to use the Dirac equation for free particles. In principle, 4'{pi,P2) is a solution of a wave 
equation which should be cast into the form of a mass operator equation by making some 
change of variable. The problem and the related constraints were considered for scalar 
particles, see ref. |27] for an example. When considering a standard meson-exchange 
interaction like the gluon one in the present case, the non-zero spin of constituents has 
to be accounted for, raising further problems. To illustrate them, it is instructive to look 
at the action of the spin part of this interaction together with that one for the pion wave 
function: 

y i^{Pl,P2) oc J2 HPiWu{p[) V{p2h^,v{p2) U{p[)-f5v{p2)- (35) 

spins 

Using a Fierz transformation, the spin part of the interaction can also be written: 

U{pih^u{p[) V{p'2)l^v{p2) = U{pi)v{p2) V{p2)u{p[) + uipi)i-f5vip2) V^p'^Yl^U^p'^) 

-^U{pi)-i^v{p2) V{p'^)-i^u{p'-^) - -u{pi)Yl5v{p2) V{p'2hf,l5u{p'l)^ (36) 

where one recognizes a term with the same structure as the input wave function (^75^) and 
another one with a pseudo- vector character {u'^^'~^^v). When performing the summation 
over spins in eq. ( !35l) . it is found that these two terms contribute, with the result: 

^ U{pi)-i^u{p\) V{p'2hf,v{p2) U{p'lh5v{p'2) = 
spins 

u{pih5v{p2) 5 — - - o npih i5v{p2) -■ 37 

2 niq 
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The first term at the r.h.s., where the dependence on variables relative to the initial and 
final states factorizes, does not provide any difficulty in reducing a wave equation to a mass 
operator. The second one does however. Writing (p'^^+pg)^ = {Pi+P2—pi—P2)fj,+ {pi+P2)fi 
and using the Dirac equation, the above equation also reads: 

u{Pih^u{p[) v{p'2)-ff,v{p2) u{p[)-f5v{p'2) = 



spins 



Pi-P'2 -r \ I \ (Pi + P'2 - Pi - P2] 



u{PihMP2) ^ - u{Pi)Y7MP2) ' (38) 

2 ruq 

where the last term is seen to depend on the surface on which physics is described (through 
^^), quite in the same way as a meson propagator would generally. The problem is well 
known in relativistic quantum mechanics. It points to the constraints that the interaction 
entering the mass operator has to fulfill. These ones implicitly account for higher-order 
corrections in the interaction, which cancel the above undesirable contributions. This is 
made possible by the fact that these contributions have an off-energy shell character (see 
eq. (jlj)). Discarding these ones, the relevant part of the gluon exchange appropriate for 
describing the pion in the RQM framework may thus be written in the following factorized 
form: 



2 U{Pl)il5U{p2) V{P2)h5V{p'l) "^g (4 + 4') \ .ggS 



The last factor has been introduced to represent as much as possible the effect of the factor 
Pi ' P'2 ^fc' + ^'^) 61- f|38|) . which originates from the original vector- vector nature of 
the interaction (7^^ 7^), while keeping the symmetry between k and k' . It differs from the 
factor (2 — mq/e|)^/^(2 — m^/e|/)^/^ introduced by Godfrey and Isgur [51] by off-energy 
shell corrections proportional to m^, which are in any case part of the uncertainty on 
the derivation of the mass operator. Due to the presence of square root factors however, 
the last choice may introduce further complication in the derivation of two-body currents. 
Whatever the choice, we notice that a pseudo-scalar pseudo-scalar type interaction (75 75) 
reproducing the non-relativistic limit of the one-gluon exchange would provide a different 
result (a factor e^/ instead of e^/ + fc'^ in the present case). Thus, the choice made here is 
consistent with the symmetry character of the interaction V ^ but could miss some off-shell 
corrections. Gathering all the relevant factors, the wave equation reads: 

(M -4e,) u{j>,)^,vij>2) - Mk) = - y J^2 24 el ) 

X—=A=u{pi)i-i^v{p2) - ^ V{p'2)i-i^u{p\) U{p[)-f5v{p'2) — Mk') 
\/ ^k'y ^k' spins 

dk' 4 ^2 I ^2 _ m2 (el + el).el 



^ u{Plh5V{P2) (2 \ 2 2 )— Mk 

fthi ^ 2 ei ei, ' mi tu 



(2 7r)3 (k-k')2 ^ 2e|e|, ^ e^. 

(40) 

which, after some simplification, is found to give eq. ([1]) (apart from the confinement part 
not considered here). 
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B High momentum behavior of solutions 



The behavior of wave functions at high momenta, which in principle is relevant for cal- 
culating form factors in the asymptotic regime, is determined by the most singular part 
of the interaction. Examination of eq. ([1]) indicates that the problem relevant to the 
determination of this behavior amounts to look for solutions of the following equation: 

(p-^)#.(r)=0. (41) 

r 

where iplr) is the Fourier transform of the momentum wave function (f){k) = y/Ck 0o(^) 
and a = Sa^/S. Though the problem was not always emphasized, it is part of earlier 
works [55l [5^ IT6]. Solutions can be obtained from the following relation: 

P{r^-') = ^^^^^^^ (r-^) , (42) 
r 

from which we infer that solutions in momentum space are given by: 

^P{p) oc , (43) 

with a = cotg{^). Particular cases are the following ones: 

a = e, fi = 1 : (p ) r ^ = 0, ipyp) oc p , 

TT r 

a = ^, /i = 1/2 : {p-^) r-^l^ = 0, ^(p) oc p-^-^ , 
2 2r 

2 2 ~ 

a = — , /i = : (p ) r^^ = 0, ipip) P"^; (limit case) . (44) 

71 Tcr 

The first case corresponds to a perturbative one which, in our case, implies the be- 
havior (j){k) oc k~^, (j)o{k) oc When the strength of the interaction increases, 
non-perturbative effects generally given by log terms can be large enough to change the 
asymptotic power, by half a unit for a = | and one unit for a = |. As has been shown 
[36] , the last case corresponds to the occurrence of a critical regime in the QED case. 
The relevance of this result for QCD and the spontaneous breaking of the chiral symme- 
try is not quite clear but we notice that the range for the strong coupling, as, usually 
referred to is of the order or even exceeds the critical value of - ~ 0.64 (a = 1.33 for 
as = 0.5). Possible problems related to the onset of a critical regime, that are generally 
ignored, appear here due to the quadratic form of the mass equation and the extra factor 
(2 — (m^ {el + e|/))/(2e| e|/)^, which separately enhance the strength of the force by a 
factor 2 at large momentum. They could be alleviated by the decreasing of as with the 
momentum transfer, or the effect of retardation effects. In practice, we will choose a value 
of as that corresponds to a = 0.5 at most {i.e. as = 3/16). 



C Two-body currents and asymptotic behavior 

We consider in this appendix the derivation of two-body currents that are related to the 
one-gluon exchange and could contribute to the pion charge form factor in the asymptotic 
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domain, of special interest in this work. This is done for the general case where physics is 
formulated on an arbitrary hyperplane of orientation and, for the following discussion, 
we introduce "positive" and "negative" energies as defined with respect to this orientation. 
Particular applications concern both the instant-form and the front-form approaches. 
Essentially, the derivation amounts to include contributions that are not accounted for in 
the single-particle contribution. This involves contributions with intermediate particles 
carrying a "negative" energy, which are not part of a RQM formalism. Contributions 
with a "positive" energy also occur but, in this case, they involve off-energy shell effects 
that the formalism ignores. These ones are not necessarily independent of the previous 
ones however. What corresponds to "negative" -energy components in a given approach 
often appears as an off-energy shell in another one. 

Our starting point is the expression of the charge form factor given by eq. fll3p . In 
order to derive the two-body currents we look for, the relation with the Born amplitude 
shown in fig. [2] has to be established. This is done by demanding that they coincide for 
the contribution with a "positive" energy intermediate quark when higher orders in the 
interaction are discarded. We notice that the above procedure fixes the strength of the 
Born amplitude entering the calculation providing a check on the ability of the underlying 
formalism to account for it. The contribution to the form factor thus obtained reads: 

dpi2 f dpf2 ^■{Pil+Pi2) flTLq 2n.\\ 



AFi(Q^) 



(2 7r)3 J (2 7r) 



9' I 



2ei2 ^-Pii 

^■{Pfl+Pf2 



^■{Pil+Pi2+Pfl+Pf2) 



(ml-pf^ f^'^-{Pi2-Pf2y 2e/2 ^-Pfi 



^{kj)) + (term : i ^ /) , (45) 



where p*j^ 



Pi2, while is given by: 



2 



xTr 



m, 



-1-Pi2 . mq+-f-pii rug+^f-pf 

^75 — 7 

2 m 2 Trig 2 irig 



r 



mg+j-pfi 
2 m„ 



«75 



-7-P/2 



2m„ 



7v 



(m^+p/i-pja) 



PiVPi2 Pf -PiVPf Pi2+Pi2-Pf Pil + mg {Pf+Pii+Pi2) 



{P*fPfl+Pf2) {Pil+Pi2) + {Pil+Pi2 ■ Pfl+Pf2) pY - {P*fPil+Pi2) (P/l+P/2) 
+ {Pi2-Pfl+Pf2) Pil - {Pil ■P/I+P/2) Pi2+PirPi2 (P/l+P/2) +"^^ (P/l+P/2) 

"(46) 

In order to make the relation with the expression of the one-body current contribution to 
the form factor, eq. ( |T3i) . the following points are noticed. The quark propagator can be 
written as the sum of two terms corresponding to "positive" and "negative" energies: 



7-p* + TUg 



Tjlg — P* 



^^il-p^ + mg) ^ {l-p +mg) 



with 



2_.*2 2i-p+ {i-p^ -i-p*) 2i-p- {i-p- -i-p*) 



(47) 



P 



'^■p* + ^-p" 
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Of course, contributions corresponding to a "negative" energy or to the part of the inter- 
action with a pseudo- vector character, which cannot be accounted for in a RQM formahsm 
and give rise to the two-body currents we looked for, are discarded at this point. Now, 
some replacements, which are justified by the lowest-order interaction terms of interest 
here, have to be made: 

9 3 3eff 3 



(49) 



2 (m^+P/rP/2) ^ iPi2+pJ)-{Pfi+Pf2) ^ fi/ _ +4/ 

m2 m'^+P~f-Pi2 ^ 

Apart from the first term at the r.h.s. of the last expression, which will be included 
in two-body currents associated with a "positive" energy intermediate state, the above 
replacements allow one to consider the integration over pf2 in eq. (H5|) . After adding the 
confining interaction, which has no direct role on the asymptotic behavior, one can use 
eq. ([T] ) and then recovers the wave function (j){kf). The identity with the expression of 
the form factor, Fi^\Q'^), eq. f|T3l) . is finally obtained by making the change of notations: 
Pi2 P, Pil Pt, Pf Pf, kf kf. 

The expression of two-body currents, which correspond to the terms discarded above, 
can now be obtained: 

(2)/ri2N _ f dpi2 f dpf2 i-{pii+Pi2) (rriq ~ \ ^(P/l+P/2) (ruq 



^ ^ i (2 7r)3 7 (2 7r)^ 2e,2 i■p^x ^e,, 2e/2 ^P/i ^e^, 



where: 



/i^ - {Pi2-Pf2y 



{ips + + Ips + Ipv) + (term : t ^ f) , (50) 



^■iPil+Pfl+Pi2+Pf2) {^■pf + ^-{Pf-p,2)) 

,^2 

X 2 



m;+Pfi-pf2 



nig 



Pii-Pi2^-Pf -Pti-Pf i-Pi2 + Pi2-Pf ^■Pii + 'ml^-{pj+pa+Pi2) 

e 



i-{Pil+Pfl+Pi2+Pf2) 2i-pf{i-pf + i-{Pf-p,2)) 
X Pj-{Pfl+Pf2) ^■{Pil+Pi2) + {Pil+Pi2)-{Pfl+Pf2) ^-Pf -pJ-{Pil+Pi2) ^(P/l+P/2) 

+Pi2-{Pfl+Pf2) ^■Pii-Pir{Pfi+Pf2) ^■Pi2 + {ml+pii-pi2) ^(p/l+P/2) 

+ = ^.(p^^ 2i.pf i^.pf - e-(P/-p.2)) 

{Pi2+p})-iPfl+Pf2) ( + 2./-+^ ^ 

^--:+ \PivPi2 i-Pt -PivPt ^■Pt2+Pi2-Pf ^■Pii + m f-{pj+pa+pi2 

m^+pJ-Pi2 \ J J J 'J J 
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-(P/ •(P/l+P/2) ^■(Pil+Pi2) + (Pil+Pi2)-(Pfl+Pf2) ^■P'f -Pp(Pil+Pi2) C-(P/l+P/2) 

+Pi2-{Pfi+Pf2) i-Pii-Pii-{Pfi+Pf2) i-Pi2 + {ml+pn-pi2) C-(P/i+P/2 



^■{Pil + Pfl + Pi2 + P/2) i^-Pf -^-{Pf -Pi2) 
1 (I 



+ 



X 



'^ ml+pfvpf2 (m^+ Pi2-P/)) + {ml+pfrPf2Y 



m„ 



m'^+ppPi2 



x{pirPi2 ^-P^ -Pn-Pj ^■Pi2+Pi2-P'f C-Pil +TOq^-(P/+Pil+Pi2)) 



(51) 



Examination of the denominator of I~ shows the presence of a factor, ^-pj + ^■{Pf—Pi2), 
which has no support in time-ordered diagrams, as can be seen in schematic models 
with scalar particles. It is cancelled by another term with the result to be replaced by 
Pf + Pii + ^'{Pf~Pi), which amounts to taking into account off-energy shell corrections. 
This has the advantage of removing a contribution that has a singular behavior in the 
limit of a massless pion (instant-form case). Thus, a more elaborate and less sensitive 
expression for I~g and which we will use, is the following: 

e 



ps 



x2 



^■{Pil+Pfl+Pi2+Pf2) '^i-Pf {^■{Pf+Pil)+^-{Pf-Pi)) 



m;+Pfi-Pf2 



^pv 



Pil-Pi2^-Pf -Pil-Pf i-Pi2+Pi2-Pf i-Pil+mli-{pf+Pii+Pi2) 

e 



X 



i-{Pii+Pn+Pi2 + Pf2) 2^-pf {^.{pf + p,^) + ^.{Pf-p,)) 

Pf-(Pfl+Pf2) ^■(Pil+Pi2) + {Pil+Pi2)-{Pfl+Pf2) ^-Pj -p]-(Pil+Pi2) C-(P/l+P/2) 



+Pi2-{Pfl+Pf2) C-Pil-Pil-iPfl+Pf2) C-Pi2 + {ml +Pil-Pi2) ^(p/l+P/2) 



(52) 



In the instant-form approach, and for the Breit frame, the above contribution does not 
show any explicit dependence on the pion mass. In the front-form approach, it vanishes 
due to the presence of the factor 

The expression for 1^^ + can also be transformed. The factor at the denominator, 
C 'Pf ~ C '{Pf~Pi2)-i can be replaced by ^- {pf + Pi2 — Pfi — P/2), which again amounts 
to neglect off-energy shell corrections. This term may be more singular but it turns 
out that the numerator vanishes at the same time. Some simplification can therefore be 
made. Moreover, the whole quantity goes to zero with the momentum transfer, allowing 
to factorize the term p'f—pa- The resulting expression we will use is thus: 

■nfi I n 



ps ' pv 



X 



^■{Pii + Pfl + P-a + Pf2) 2^-pf {ml+pf ■pi2) 
C-P~f {{P~f-Pii)-C {ml+Pii-Pi2) - {Pii+Pi2)-C {P'f-Pii)-Pi2) 
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-(P/+Pi2)-C Pi2-i {P'f-Pil)-Pil 

+ {{p'f-Pn)-{Pii+Pi2) - {p'f-PiiX Pi2-i) {ml+p'^-pi2) 

e{p}-P^l)< 

-l 

C • {Pil + Pfl + Pi2 + P/2) '^i-Pfi- {Pil + + 2 Pi2) C ■ {Pi2 +Pf+ Pfl + P/2) 

^ \L ml+Pfi-Pf2 _ {ml+Pi2-pp) + {ml+Pfi-pf2) \ 
[\ rnl rnl+pppi2 J 

x{pii-Pi2 C-P'f -Pii-pj ^■Pi2+Pi2-P'f ^■Pii + ml^-(pj+pa+pi2)) . 

(53) 

It can be checked that the above expression vanishes for both the instant form and the 
Breit frame (Pj = —Pi). It does not however away from this configuration. In the 
configuration Ei = Ef together with the limit |Pj + Pj| 00, the expression is found to 
be identical to the usual front-form one (g+ = 0). This one could be obtained from the 
above one using standard changes of variable. 

The derivation of two-body currents considered in this appendix could be easily ex- 
tended to the Lorentz-scalar form factor. As this one already evidences the expected 
asymptotic power law, except perhaps for a numerical factor, the need for considering 
two-body currents is not as strong as for the charge form factor. Moreover, a perturba- 
tive calculation may reveal to provide little effect in the domain where it does for the 
other form factor. 
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